Quantum Noise, Scaling and Domain Formation in a Spinor BEC 
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In this paper we discuss Bose- Einstein spinor condensates for F=l atoms in the context of ^'^Rb, 
as studied experimentally by the Stamper-Kurni group. The dynamical quantum fluctuations of 
a sample that starts as a condensate of A'' atoms in a pure F = 1, mp — state are described 
in analogy to the 'two-mode squeezing' of quantum optics in terms of an su(l,l) algebra. In this 
system the initial mp = condensate acts as a source ('pump') for the creation pairs of mp =1,-1 
atoms. We show that even though the system as a whole is described by a pure state with zero 
entropy, the reduced density matrix for the mp = +1 degree of freedom, obtained by tracing out 
the mp = -1,0 degrees of freedom, corresponds to a thermal state. Furthermore, these quantum 
fluctuations of the initial dynamics of the system provide the seeds for the formation of domains of 
ferromagnetically aligned spins. 

PACS numbers; 03.75.Mn, 03.75.Lm, 03.75.Kk 



I. 



INTRODUCTION 



Spinor condensates were first realized^i^ in 1998. Such 
condensates are very rich in the underlying physics 
and have been the subject of numerous studies, from 
mean field theory and ground state considerations^!^ 
to experimental observation of dynamics, rotation ef- 
fects, domain structures^i^ '^^i^^'^^ , theoretical dynam- 
ics of mixing^-ii^ii^; existence of quantum vortices and 
topological structurei^iiiiiSiii^; dipolar effects^S and very 
many more, currently making this one of the most ac- 
tive fields at the interface between condensed matter and 
atomic physics. 

In this paper we investigate the formation of domains 
in terms of a Bose-Einstein spinor condensate of alkali 
atoms, ®^ Rb, in the hyperfine multiplet F — 1, which has 
been shown to have ferromagnetic interaction amongst 
the constituent atoms^i^i^. A recent experiment carried 
out by the Stamper-Kurn group at Berkeley-'^, showed 
the formation of domains in a ®''Rb gas starting with 
a polar non-ferromagnetic initial state. The experimen- 
tal sample formed random domains of varying transverse 
magnetization, and this paper is based on considerations 
regarding this experiment and the dynamical seeding of 
the observed domains. The analysis of this experiment 
has been undertaken by Lamacraft^^ who derived a form 
for the dynamics of the particles in different states. The 
dynamics of the associated quantum phase transition and 
scaling properties of the system magnetization has been 
discussed by Damski et al.'^^ for a one-dimensional sys- 
tem treated in a mean-field approach, with extensions 
to higher dimensions. The recent work of Saito et ali^ 
further investigated the dynamic formation of these do- 
mains in terms of quantum noise, including a simulation 
of two-dimensional domains using a combination of quan- 
tum dynamics and Gross-Pitaevskii evolution - this work 
was completed concurrently with our own investigation 
with essentially the same conclusions^'^. In this paper 
we discuss an alternate, yet equivalent description of the 



experimental work, which elucidates the nature of the 
quantum statistics. After introducing the system and the 
relevant experiment in Sees. HlHIl we proceed in Sec. IIIII 
to supplement previous results on the single-mode prob- 
lem by providing a connection to the ideas of quantum 
squeezing and quantum noise, as well as some exact solu- 
tions to the dynamical single-mode problem. In Sec. IIVI 
we also follow a multi-mode approach, in terms of which 
the time evolution, quantum noise and statistics of the 
domain seeding may be discussed. This is presented in 
connection to an su(l,l) algebra inherent in the effective 
Hamiltonian derived for the system, making it possible 
to obtain probability distributions for the fluctuations in 
the system. As we will see the thermal nature of these 
distributions is distinctive and might be verifiable exper- 
imentally. 



A. Effective theory of spinor condensates 

The use of dipolar optical traps allows for Bose- 
Einstein condensation of alkali atoms in which the spin 
degree of freedom is still active - the traps do not pref- 
erentially select one of the spin states with a given mp 
(as opposed to magnetic traps that favor the weak-field 
seeking states) . Thus for a gas of spin F there are 2F + 1 
degrees of freedom deriving from the hyperfine spin. The 
total spin for both ^^Rb and ^■^Na is the sum of the nu- 
clear spin, / = 3/2, and the electronic spin 5* ~ 1/2, 
leading to a total spin F = 3/2 ± 1/2. For the F = 1 
manifold, which is usually probed in optical trap exper- 
iments, this means that the atoms may be described by 
a spinor, 



Mr) 
Mr) 
M{r) 



(1) 



with each component, tpi{r), corresponding to the wave 
function of a species in the mp = {—1,0, 1} state. At 



2 



zero magnetic field the system's spin rotational invari- 
ance is manifest. The standard Hamiltonian for the low 
energy dynamics of an i<" = 1 dilute atomic gas was de- 
veloped by Hc^ and Ohmi and Machida^. This assumes 
that only two-body collisions are important and that the 
atoms do not interact otherwise with each other and that 
the system is rotationally invariant, which means that the 
interactions can only depend on the total spin T of the 
colliding atoms and not on ^ ^24^25,26^27 ^ complete 
effective Hamiltonian for the F = 1 dilute cold gases in 
second quantized form is written as, 



H = J d^r|v'l(r) 



2m 



+V'l(r)l/yV',(r) + i 



co-0l(r)V'!(r)'0*(r)-0j(r) 



+C2 (vl(r)F,,V,(r)) • (v4(r)Ffei^Kr) 



, (2) 



where Vij is the trapping potential, the indices refer to 
the hyperfine species, i,j £ {—1,0, 1}, and we are using 
the Einstein summation convention of summing over re- 
peated indices. F = {Fr^, Fy, F^} is a vector of spin-one 
matrices. We have also identified the couplings 



Co 



(2a2 + ao), 



3m 

C2 = -z — (02 - ao), 
6m 



(3) 
(4) 



which depend on only two parameters, the scattering 
lengths ao and 02. The field operators obey the com- 
mutation relations. 



(5) 



with all other commutators being zero. The trapping po- 
tential, Vij is the result of a combination of magnetic and 
optical potentials and may be taken as diagonal, V{r), 



assuming that all magnetic fields are in the 'z-direction'. 
We should keep in mind that in writing the Hamiltonian 
in terms of projecting onto total spin f subspaces we are 
assuming that any Zeeman shifts are smaller than the hy- 
perfine splitting so that J- is still a good quantum num- 
ber. The mean field analysis of this Hamiltonian^'^'^ii^ 
indicates that the ground state for C2 < 0, as is the case 
for ^^Rb, is ferromagnetic. 
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(6) 



which is written this way to emphasize that we have a 
degenerate set of ground states which are related to each 
other by a gauge transformation, exp[j0], and an arbi- 
trary rotation via 

C/(/3i,/32,/?3) = exp[-iF,/3i]exp[-zFj,/?2]exp[-iF,/33], 

(7) 

where {Pi, (32,^3} are the Euler angles^^. The applica- 
tion of a magnetic field would lead to the average spin 
pointing in the direction of the applied field and the se- 
lection of one of the degenerate ground states. In the ab- 
sence of a field we have spontaneous symmetry breaking, 
where the system selects at random one of the possible 
directions for its spins. 

We may write the field operators for the case of a ho- 
mogeneous system as a plane wave expansion: 



ipi{r) = a,k exp k • r] . 



(8) 



This gives: 



Hi = 



X! I (co + C2)(alkialk2aik3aik4 + a-iki«-ik2«-ik3a-ik4) + coaJ)kiaSk2«ok3aok4 



ki+k2=k3+k4 



+2(co + C2)(aik^aJk^aik3aok4 + a-ikiaok2^-ik3aok4) + 2(co - C2)a\:^^^a}_^:^^^ai-k^a-i)^^ 
+2c2(alkia-ik2"ok3aok4 + aJkiaok2aik3a-ik4)|- 

I 



(9) 



The interaction involves self-scattering and cross- 
scattering terms among the three particle flavors, in ad- 
dition to the last term which involves the conversion of 
pairs oi mp ~ Q particles to mp = ±1 and vice versa. 
This provides for interesting interspecies dynamics that 
we will explore further in later sections. 



II. AN INTERESTING EXPERIMENT ON **^RB 



In an interesting experiment carried out by the 
Stamper-Kurn group at Berkeley^, an initial sample of 
^''Rb atoms was prepared in the mp = — 1 state and 
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trapped in a quasi-two-dimensional optical trap with os- 
cillation frequencies {ijjx,ujy,ujz} = 27r{56, 350, 4.3}s~^ 
at a longitudinal magnetic field = 2G. Subsequently, 
the atoms were converted to the mp = state by the 
application of a radio frequcncy(r.f.) field, reaching a 
peak density of n = 2.8 x lO^^cm^'^. The magnetic field 
was then quickly ramped down linearly in 5ms to about 
50mG and the gas was allowed to evolve freely in the 
trap. The presence of the original magnetic field pro- 
vides a quadratic Zeeman interaction which lifts the de- 
generacy for the transitions | = — 1 ) — > | jnp = ) 
and Imp = 0) ^ Imp = I) states, so that good conver- 
sion may be achieved from the initial mp = —1 states to 
TTiF = states. Since C2 < for ^^Rb, the interactions 
are ferromagnetic in nature, and it is thus energetically 
favorable for the spins to align with each other. So, when 
the magnetic field is ramped down and becomes effec- 
tively zero, the mp — states become dynamically un- 
stable because of the exchange interaction and convert to 
mp — ±1 pairs. The different variables in the experiment 
are summarized in Table(|T]). Within some time, Thoi^^ the 
transverse magnetization in the a^y-plane was then im- 
aged, using a novel non-destructive in situ technique, by 
detecting its Larmor precession about a guide field. The 
experiment revealed that the longitudinal magnetization 
was negligible. In contrast, the images of the transverse 
magnetization, indicated the formation of multiple ran- 
domly oriented domains of varying shapes and sizes, as 
well as more involved spin textures. The typical size, 
^exp, of the domains seen after the domain growth satu- 
rated was in the region of ~ lO/im. The growth of the 
transverse magnetization was observed to be initially ex- 
ponential, with a time constant t ~ 15(4)ms. 



If one tries to model the initial seeding of these do- 
mains using the Gross-Pitaevskii equations and starting 
from a pure | N)^^^^ state, then no evolution of trans- 
verse magnetization is observe d^^'^^ . We may think of 
the Gross-Pitaevskii equations as corresponding to non- 
linear Schrodinger equation describing the motion that 
begins sitting on top of a potential hill. Without an 
initial displacement nothing happens classically. Some 
form of noise is required to provide the instability to roll 
off down the potential hill. In a recent paper by Saito 
et alJ^ different forms of noise, such as white noise or 
noise due to the unconverted initial mp ^ —1 popula- 
tion, were tried out to reproduce the instability seen in 
the experiment by Sadler et al.^. In this paper we in- 
stead model the instability that causes the seeding of the 
domains observed experimentally in terms of quantum 
noise, and obtain an analytic form for this in the initial 
time regimes that the seeding takes place. Recent work 
by Saito et alj^ uses a similar picture. 



'^'Rb data 


ao scattering length 


101.8 as 


a2 scattering length 


100.4 as 


Mass 


87 g moP^ 


Experimental Parameters 


Number of atoms 


2.1(l)xl0'' 


Peak density, n 


2.8x10" cm"^ 


Temperature, T 


40 nK 


con/ks energy scale 


~ 100 nK 


C2n/kB energy scale 


~ 480 pK 


Trap frequencies {ujx,i^y,i^z} 


27r{56,350,4.3} s^^ 


Initial magnetic field, 


2 G 


Final magnetic field, Bz 


50 mG 


Time duration, Thou 


36-216 ms 


fi/(c2n)time scale 


~15.8 ms 


Observed time constant, r 


15(4) ms 


Typical domain size, £,exp 


~10 /im 



TABLE I: Experimental values of interest for an experiment 



III. THE SINGLE MODE HAMILTONIAN AND 
ITS DYNAMICS 

The simplest approach one can take beyond mean-field 
theory is to look at the single mode Hamiltonian. This 
has been done by various authorsi^i^S, including some 
numerical work within a classical frameworklSi^i^. In 
this section we will independently reproduce some of the 
previous results in a different approach and supplement 
them with new insights. 

We start with the effective Hamiltonian, for N spin- 
1 bosons in zero magnetic field and we assume that all 
the mp species have the same spatial wave function, 
77(r). This wave function is determined by using a Gross 
Pitaevskii equation which comes from the kinetic part 
of the Hamiltonian and the cq coupling - the symmetric 
parts of the Hamiltonian, Eq.Q. Namely, 

+ ^ + coiV|77(r)|2^ rjir) = fiovir), (10) 

where fiQ is a chemical potential enforcing the total par- 
ticle conservation. Since for ^''Rb cq 3> C2, this means 
the spatial wave functions, ?7(r), of the condensate are 
largely dependent on cq and may be taken to be the same 
for each mp species to a first aproximation. The fields 
are thus approximated by 

V'j(r) = ?/(r)ai; i/',(r)'f = vi'^)*al. (11) 

Here we have defined the operators a^, al that respec- 
tively annihilate or create a boson in the state = 
mp — i, with i G {—1,0,1} in the z-basis representa- 
tion of the spin-1 operators. These operators obey the 
commutation relation 

= '^M' (12) 

with other commutators being zero. We can further de- 
fine a number operator that counts the number of bosons 



4 



in the state i as Ni = ajoi. Using these operators, 
and taking into account the total number conservation, 
iV = iVo + iVi + 7V_ 1 , we may rewrite the Hamiltonian in 
the form of Diener et a\^: 



Ho = -( (iVi-iV_i)2 + (2A^o-l)(A^i+A^-i) 



t^2 



c is given by: 



C2 / d'r\r,{r)\' 



(13) 



(14) 



Notice here that the cq terms which would relate to den- 
sity fluctuations vanish up to a constant and thus do not 
affect the dynamics of the system and do not enter the 
single-mode Hamiltonian considered above. 



A. Dynamical considerations 



This system of linear differential equations can be solved 
exactly to give: 

ai{t) = Ai{t)ai + A^i{t)al,, (18) 
a-i{t) = A{{t)al^+ A*_^{t)ai. (19) 

where operators with no explicit time dependence corre- 
spond to operators at initial times, t = 0, and in addition 
we have defined. 



Aiit) 



cosh 



cWANa - 1 



2h 



A-,{t) ^ 



-li^^sinh 
2iNo 



cty/ANo - 1 

2h 



sinh 



cWiNo - 1 



2h 



(20) 
(21) 



We make here the Bogoliubov approximation that we 
may replace the annihilation-creation operators for the 
condensate with numbers, oq « aj « y/ Nq, up to a phase 
factor that we may neglect since we are later concerned 
only with expectation values where the phase would can- 
cel out. This approximation is valid for extremely short 



times as if our initial state |A'^)q condensate acts as 
an unlimited source for the creation of particles in the 
tof = ±1 states that does not get depleted. Also, we 
notice that in this approximation, the numbers of par- 
ticles in the two spin states mp = ±1 are equal at all 
times. This means we can ignore the first term in Eq. (jl3p 
that involes the difference {Ni — N-i). Then the Bogoli- 
ubov approximation Hamiltonian becomes quadratic in 
the fields, 

Hb^^ [{2No - l){Ni + iV_i) + 2Noia\aU + aiO-i)) • 

(15) 

We may obtain Heisenberg equations of motion in this 
regime as follows: 

ihdtai = [ai,i/o] = ^{2No - l)ai + cNoal^; (16) 



B. Number averages and variances 

It is interesting to calculate the expectation value 
of the operator Ni{t) starting with a ground state 
where A^i — 0. From the results in the pre- 
vious section and starting with the ground state 
(OU=_i(A^L^.o(OU=i ^ (0;^;0| we get the 
number expectation values, 

N^{t) = {0;N;0\N^it)\0;N;0) 



2(iVo)sinh 



ct^/4No-l 
2h 



y/ANo - 1 



(22) 



'-{2NQ^l)al^-cNoai. 



(17) 



with /i G {±1}. This indicates rapid population transfor- 
mation from A'o to equal numbers of mp = ±1 particles. 
The equality oi mp populations is guaranteed from the 
Hamiltonian and the conservation of the total number 
of particles and the total spin in the system. We may 
investigate the single-mode dynamics further by calcu- 
lating the variance. 



N,{tY-(N,{t)f ^\A,\M-i\ 

(^0 si 



4 sinh 



ct^/ANg-l 

2h 



sinh 



2h 



N, 



(1 - 4A^o)^ 



(23) 

(24) 
(25) 
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Scaled time evolution for various N 
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FIG. 1: (Color online) The single mode Hamiltonian may 
be solved exactly for different A^. The time evolution of the 
solutions indicates scaling of t as the curves collapse 

into a single line for short times. 



Unsealed time evolution for various N 




FIG. 2: (Color online) The exact solutions for the single mode 
Hamiltonian exhibit a linear time regime, with the curves for 
different A'^ collapsing into a single curve without scaling. The 
populations oscillate and ultimately complete many cycles of 
non-linear evolution followed by linear evolution. 



This is exactly the variance of a thermal stated and it 
actually corresponds to a Bose-Einstein distribution - we 
will see this explicitly in the context of our multiple mode 
treatment in Sec. (|IV C[) . The single mode condensate 
starts out with a source of N particles in the mp — 
pure state that are transmuted into equal numbers of 
mp = ±1. Each of the resulting species displays a super- 
Poissonian variance, having a temperature and entropy 
associated with it. We should keep in mind that this is 
done in the Bogoliubov approximation and is in need of 
further exploration. Such considerations are undertaken 
in the next section for exact single mode solutions for 
various numbers of particles, as well as in Sec. pvp . which 
explores a multi-mode approximation. 



where < i,j < N/2. We may think of this as a hopping 
Hamiltonian, allowing for movement between adjacent 
states that differ by one pair. This is essentially a one- 
dimensional problem, resulting in a tridiagonal matrix 
that may be solved numerically to obtain the eigenvalues 
and time evolution of the system. The time evolution for 
various N is shown in the two figures, Fig.(IT|) and Fig. ([2]). 

As we can see from the figures, we may identify a crit- 
ical time tc ~ which separates two time regimes, in 

agreement with the results of Law et a]M. We also no- 
tice the additional feature of scaling with for short 
time scales, where the different evolutions for different 
total number of particles collapse into a single curve as 



C. Exact single mode solutions for different 
numbers of particles 

We now turn our attention to solving the effective sin- 
gle mode Hamiltonian, Eq. exactly for different ini- 
tial numbers of particles N. Given our starting state 
I -^)mF=0' O'^ly pairs of particles may be created or de- 
stroyed, and this makes for a smaller subspace under our 
consideration. If we take our initial state to be the 'num- 
ber of pairs' vacuum state, for A'^ particles we may only 
have Y + 1 number of pair-states. We may write the ma- 
trix components of the Hamiltonian that connects pair- 
states with i and j pairs respectively as: 



H, 



_ -^2j(2(7V-2z)-1)<5,,, 



+2^{N-2j){N-2j-l)S,^j+^ 
+2 V(iV - 2j + 1){N - 2j + 2)5,,,_i }, (26) 



- 1 

O.B 
0.6 
0.4 
0.2 



Exact Vs. approximate solutions 



/ 



Bogoliubov 
N =204 




FIG. 3: The exact solution for A'^ = 204 is compared to the 
Bogoliubov approximation. As is seen from the graph the 
two are in good agreement at the very early initial times, 
in support of our no-depletion approximation for the single 
mode Hamiltonian. 
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Time-scaled deviation ratio for various N 



Deviation ratio for various 



ratio 




FIG. 4: (Color online) The ratio ^ N^,{t){N^{t) + 1) to the 
deviation cr is a measure of the correlation of the state variance 
to that of a thermal state. Our Bogoliubov approximation is 
shown to obey a thermal variance and to match the exact 
results only for very short times. The various curves can be 
seen to collapse again into a single curve showing that the 
plotted ratio is independent of A''. 



seen in Fig.([T]). There are two time regimes in the sys- 
tem: First a regime within which time scales as when 
the dynamics is dominated by the non-linear part of the 
Hamiltonian shown in Fig.((T]). Then we see a subse- 
quent metastable time regime, when the linear part of 
the Hamiltonian takes over the time evolution and the 
time does not scale with the number of particles, Fig. ([2]). 
The time regimes alternate for many cycles if the system 
is allowed to evolve over time. The starting condensate 
of N particles in the mp — state is initially rapidly 
depleted to create pairs of mi? = ±1 particles. When 

the metastable time regime is reached, at i^^^ = 1, the 
population ratios are 



7V_ 



N 



Ni No 
N ' N 



0.25 : 0.25 : 0.5. 



(27) 



This metastable regime is followed by oscillations in the 
three different populations and eventually the complete 
repopulation of the mp = state, reaching again the 
initial state and the cycle repeating itself when — 2tt. 

Furthermore, we may compare our Bogoliubov trans- 
formation for short times to the exact solution for a 
given N. As shown in Fig. ([3]) the two are in agree- 
ment only for very short times. Assuming the scaling 
of time with this implies a validity time dictated 
by the depletion, where the percentage, kq, of parti- 
cles remaining in the | N)^^^^ state is on the order of 



No 



0.9. We are in- 



terested in the seeding of the domains at very early times 
in the time-evolution and hence the Bogoliubov approx- 
imation looks indeed like a good starting point. 

Additionally, the statistical variance and deviation 
were investigated for various functional forms of the num- 



ratio 




FIG. 5: (Color online) The ratio \/ N^(t){N^{t) + 1) to the 
deviation a is here shown for long times. We still observe 
that our Bogoliubov approximation always obeys a thermal 
variance. Furthermore, the plot indicates that the plotted 
ratio is independent of indicating that the actual variance 
is some functional of a thermal variance. 



ber expectation, A'^ (i). It was found that the curves for 
different number of particles collapse into single curves 
when we scale the variance with N^{t){Nfj,{t) + 1). Plots 
for different time scalings corresponding to the two time 
regimes identified above are shown in Fig. ([4]) and Fig. ([5]). 

In the plots the ratio of ^JW^^(W^^'+l) to the de- 
viation <T, reveals scaling and collapse of the plotted 
curves showing that the deviation must be intimately 
related to a thermal distribution. The results of the 
Bogoliubov approximation are included for reference. 
As shown, the Bogoliubov approximation obeys a ther- 
mal variance through all time regimes, displaying super- 
Poissonian statistics consistent with two-mode squeezing 



picture" 



33.34.35.36 



I, which will be discussed further in terms 



of the multimode approximation in the next section. 



IV. MULTIMODE BOGOLIUBOV 
APPROXIMATION 

We now would like to extend our discussion to the case 
of many spatial modes. The picture we have in mind is 
that mp = ±1 pairs are created from the mp = source 
and then separate to begin seeding of the domains. The 
interplay of the various wavevectors k will then be im- 
portant in ascertaining the real space distribution and 
behavior of these domain seeds. We first look at an ef- 
fective Hamiltonian description. 



A. The effective multimode Hamiltonian 

We begin with Eq.® and make the approximations 
that: we have no mp = atoms in a non-condensed 
state; self scattering and cross scattering terms outside 
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of the spinor mp- = condensate maybe be ignored; we 
may take our Bogoliubov approximation once more for 
the niF — operators, ao 
Hamihonian becomes, 



N. Then the effective 



HiB = 



2V 



(cq + C2) 

V 



(alk«ik + alikfl-ik) 



+ ^(«ikaii-k + aikfl-i-k) 



(28) 



As a simphfication we will be assuming we can drop the 
Co terms in some short-time 'no-depletion' approximation 
since the variation of the cq terms with respect to A^o ap- 
proaches zero in this approximation - as we recall the 
terms cancel out identically in the single-mode approxi- 
mation. The magnitude of the coupling ratio ^'--^10'* 
suggests that density fluctuation modes, as dictated by 
the Cq terms, are very stiff, so they are ignored in our cal- 
culations. We further neglect the spatial dependence of 
the trap potential. With all these simplifications made, 
we can finally write down the effective Hamiltonian we 
will use for our multimode treatment of the spinor con- 
densate. 

He = ^ s(ek + C2n)(a|kaik + alikfl-ik) 
1^ ^ 



+C2n{a{^^a^_-^_^^ + aika_i_k) } , (29) 



(30) 



where we have defined the kinetic energy, 

The effective Hamiltonian above treats the system in 
terms of pairs of the two species with mi? = ±1 having 
opposite wavevectors. The first term in the Hamiltonian 
counts the number of pairs and the second term provides 
an interaction via creation and annihilation of the said 
pairs. The Hamiltonian is quadratic in the field opera- 
tors and thus we are looking fluctuations of the system 
at a gaussian level. 

B. Time evolution of operators, noise and statistics 

For a given mode we may obtain the following equa- 
tions of motion, 

(9ciik(^) , \ / \ t / \ 

ih — — — = (ek + C2n)aii^(t) + C2na'_^_-^(t); 



-ih 



dt 

9aVkW 
dt 



(ek + C2n) a^i_k(*) + C2"aik(0- 



(31) 



As for the single mode case these may be solved exactly 
to obtain 

aik(i) = ^ik(i)aik + ^-i-k(^)ali_k; 
a-i-kW = ^lkWa^-k + ^-i-k(i)aik, (32) 



where we now have defined 
t 



Aik(t) = 



cosh 
—ic2n 



sinh 



i{c2n + ek) . 

■ smh 



Gk 



:Gk 



(33) 



with an effective 'gain' function for a given mode k, 

G2 =ek(-2c2n-ek). (34) 

Again, in these solutions operators with no explicit time 
dependence refer to the inital time, t = 0, operators. The 
paper by Saito et ali^ also arrives at the same solutions 
to the differential equations, with a similar treatment of 
ignoring the Co terms. The current work supplements 
these with single mode considerations, and with an al- 
ternative treatment that is developed below: a quantum 
noise approach where a new connection is made to the 
5u(l, 1) language of quantum squeezing. 

The two operators aik, a_i_k are closely intertwined, 
with the population of one acting as a noise source for the 
other and vice versa. This suggests also that the growth 
of different species with opposite wavevectors is a highly 
correlated process. Indeed, the correlation of species of 
opposite wavevectors is non-zero, 

™,=o ( N I a-ik'(t)aik(t) I A^>™,=o = 

n 2 



C2n(c2n + ek) . , 
- smh 



iC2n 

"gT 



Gk 



sinh 



:Gk 













cosh 




} 











^4,-k', (35) 



where we may picture that the creation of a 1 1 )]^k state 
is correlated with the creation of a | l)_i_y^ state. The 
A^-particle condensate source acts as a reservoir for the 
creation of pairs of particles moving in opposite direc- 
tions. The mp = ±1 states are populated rapidly, at 
an exponential rate, as shown by the evolution of the 
number expectation for a given wavevector: 



.=o{N\N±i±k\N' 



rriF—O 



■ sinh 



We should note here that 



limAr^k(i) 



cWt^ 



(36) 



(37) 



which agrees with our single mode considerations in the 
same limit, Eq. (P^ . 

The exponential growth of the number of particles with 
time implies our no-depletion approximation is only valid 
in the initial stages of the time evolution. Furthermore, 
the growth is dominated by certain values of the magni- 
tude of the wavevectors k as seen by the behavior of G|, 
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Dominant fc values 



Depletion of initial source 




■\/2mlc2\n 




FIG. 6: The gain parameter Gk is real for only a set of values 
of k and these are the major contributors in the exponential 
growth of the number of mp = ±1 particles. Imaginary values 
lead to oscillatory behavior instead which does not support 
the number growth of the two species. 



FIG. 8: The population of the mp = ±1 species leads to the 
depletion of the initial source of A'o atoms in the mp = 
state. Since our approximation assumes no depletion we can 
only ascertain its validity in the initial time region in the 
figure where it is only a few percent of the total available A^o 
reservoir. 



Number evolution 




FIG. 7: (Color online) For ^''Rb, the expectation number of 
particles for a given species grows over a certain range of val- 
ues of the wavevector magnitude k. Furthermore the growth 
in time is exponential and so our no-depletion approximation 
remains valid only at earlier times 



Fig.Q, and is independent of direction. The dominant 
wavector is at the maximum of G^: 



h — 



-2mcnn 



fi2 



(38) 



This is the wave vector with the greatest initial insta- 
bility. For *^Rb, C2 < and hence kmax € R- It defines 



a coherence length scale, 



Co 



2mc2n 



2.4^m, 



which in turn gives a coherence volume 



^0 



(39) 



(40) 



which for the experimental peak density, n = 2.8 x 
10^ ^"^cm^^, has associated with it = nV^g ~ 1.6 x 10* 
particles. This is large enough to support our assump- 
tion that the initial fluctuations are gaussian in nature, 
as suggested in our approximation of using a quadratic 
Hamiltonian. This allows us to explore the small quan- 
tum fluctuations that provide the starting instability for 
the evolution of the system - assuming no depletion of 
our initial particle source. We may investigate further 
by looking at the number evolution in time as a function 
of wavevector as shown in Fig. ([7]). This again shows the 
range of wavevectors over which the instability causes 
the growth of the mp = ±1 states. Additionally, we 
observe that the growth is exponential in time, and that 
initially our depletion approximation should be valid. As 
an aside, we note here that for ^'^Na Gk does not exhibit a 
real kmax since for this system C2 >0. Hence for ^^Na the 
number expectation values are oscillatory and no strong 
evolution into mp — ±1 species would take place, con- 
sistent with the 'antiferromagnetic' interactions in this 
systemSI. We may calculate the percentage of particles 
remaining in the state source of the system as 

a function of time. We have. 

No • i J 
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The sum may not be calculated in closed form, but may 
be approximated by integrating the wavevector over the 
"instability" range. The results are shown in Fig.®, 
again indicating that at early times the depletion is only 
a few percent and our approximations must still be valid. 
For the characteristic time t^^^^ we can substitute for 
the Berkeley experiment parameters to obtain each time 
unit as ^ 16ms. Then up to say 4 x 16ms in Fig.® 
corresponds to the first two initial images from the ex- 
periment, as shown in Fig. (2) in the paper by Sadler et 
al.^ . The experimental results are consistent with our no- 
depletion approximation for the initial times, when the 
domains are still forming and the transverse magnetiza- 
tion has not yet saturated as seen from the brightness of 
the figures. 



1. Thermal Statistics 

Having calculated the number expectation for a given 
species and mode it is natural to consider its variance. It 
turns out that we may express our result simply, in terms 
of the number operator expectation value in the mp = 
starting state evolving in time. 



+ c^v? sinh^ 



Gk])sinh2 [iGk 



Gi 



(42) 



This is exactly the same form as for a thermal state. 
This may be slightly counterintuitive because we start 
out with a pure state at zero temperature, with no en- 
tropy assocciated with it and end up with essentially a 
thermal distribution (albeit with a different temperature 
for each spatial k mode) . If we were to imagine that we 
were aware of the existence of only one of the two species, 
and carried a measurement of only that species for a given 
wavevector, we would associate a temperature to it based 
on its variance. Equivalently, we may think of this as 
what we would get if say we traced out the mp—-l atoms: 
the m_F=l atoms display an effective finite entropy and 
a thermal distribution. These super-Poissonian statis- 
tics show atom bunching and bear similarities to pro- 
cesses in as diverse fields as quantum optic o'^^'^'^i'^^'^^ , 
topological defects and particle production in cosmology 
and gravitatio n^^i^^i^^i'*°i'^^i"^^i^'^ , within the common lan- 
guage of squeezing and the su(l,l) algebra with which we 



will recast our results in the following section. Further- 
more, we will see how the thermal variance corresponds 
to a Bose-Einstein distribution using a reduced density 
matrix computation. 



C. The su(l,l) underlying algebra and evolution 

We notice that our Hamiltonian Hp has an underlying 
structure in terms of pairs of different species and oppo- 
site wavevectors. Namely, we define operators that anni- 
hilate and create pairs at a given wavevector k, and 
respectively, as well as an operator counting pairs, 



Ky = aika_i_k, 
K+ = atk«-i-k' 

= 2 ("ik^^ik + a-i-kfl^i^k 



(43) 
(44) 

(45) 



The commutators of these operators turn out to be: 

[K,K^] = ±Kt (46) 



[^k,^k] - -2i^k, 



(47) 



providing a realization of the generators of an su(l,l) 
algebrai^*^. In terms of these we may rewrite our effec- 
tive Hamiltonian as, 



He - ^ i (ek + C2n) {2K^ - 1) + nc2{K+ + if^ 



(48) 

We observe that the operators create and annihi- 
late pairs of atoms of different mp species, with oppo- 
site wavevectors. Then simply counts the number 
of these pairs in the system. We have thus arrived at 
the su(l,l) algebra description that is associated with 
quantum squeezing. This is a multimode realization of 
the algebra, one realization per pair of atoms of dif- 
ferent species and oppositely directed wavevectors. As 
we have already seen in the previous section the associ- 
ated quantum states are characterized by thermal, super- 
Poissonian fluctuations, as was already discovered in our 
single-mode considerations. 

Having obtained the form of the Hamiltonian we may 
now use it for unitary time evolution to arrive after time 
t to some final state | / ) , 



/) = exp 



-iHst' 



\n: 



mp-—0 



[ -it ( 

Yl exp[i0k] exp — <^ (ek + czn) 2K^ + nc2{K+ + 



rriF—O ' 



(49) 



(50) 
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where, 

t 



k = T (ek + C2n) . (51) 
n 



We want to disentangle the exponential as: 

exp[0{a+A'+ + aik^^ + a^K^}] - exp[0+(0)7^+] eM4>l{0)Kl] exp[0^(0)i^k ], (52) 



where 9 is an auxiliary parameter which is set to one at 
the end and we have defined: 



-it 
n 



(53) 
(54) 



Matching the coefficients of the generators on both sides 
would give^^i^: 



^k - 0k exp[- 



ill 



ill 



(55) 
(56) 
(57) 



The last two equations may be substituted into the first 
to obtain a Ricatti equation^^, 



^+al 



0, 



(58) 



which gives the solutions: 

sinh(rk6') 



Km 
0k (^^) 



r 



k cosh(rk6') - 
2 ln[cosh(rk6 



a° sinh(rk6l) 



(59) 



2rk 



^sinh(rk0)], (60) 



sinh(ri< 



cosh(ri. 



a° sinh(rk9) 
2IV 



where 



-4«k 



— r2 

;^2^k- 



Our final state becomes: 



(61) 



(62) 
(63) 





\K 1 


|exp 





;)"i^>iki«)-i- 



®ki/k(t); 



(64) 



This is an exact solution for the dynamics of our effec- 
tive Hamiltonian, He, in terms of number states for a 
given wavevector and particle species. We can either use 
the time-evolved state or simply the time-evolution of the 
creation annihilation operators to investigate the statis- 
tical behavior of the system and arrive at the familiar 
expression for two-mode squeezing 

(A^pk(i)') - (iV^kW') - {N^^{t)) {{N^^it)) + 1) , (65) 

once again indicating that the variance for a given mode 



is thermal in nature. 

Let us investigate the thermal variance a bit further. 
Say we look at the density operator pk for a given 
wavevector, which for a pure state corresponds to 

Pk = |/k(t))(/k(i)l- (66) 

We now perform a partial trace, to obtain the density 
matrix for the mp = 1 particles after tracing out the 
mp = —\ degrees of freedom for a given wavevector k, 



p+ik = Trinp=_i/5k 

n2cisinh[^]' 



E 



cosh [^] ' + (^i^ sinh m'^\ Gl (cosh ' + sinh ' 



I '/Ik \ Mlk 



1 ivik(o Y 

l + iVik(t) V Vi + ^ikW / 



/ik(*lik- (67) 
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An identical calculation but instead tracing out the mp = I states gives, 



P-ik 



r 



The reduced density operators, p±ik, have exactly 
the form of the density matrix for a Bose-Einstein 
distributional^. This is remarkable because our pair 
state I /(i) ) is a pure state that has no entropy associated 
with it. On the other hand the reduced density operators, 
P±ik, correspond to thermal Bose-Einstein distributions, 
that have an associated entropy and noise that results 
from the perfect correlations between the mp — ±1 par- 
ticles in each pair state. For a given wavevector, the 
mp = -l-l particles act as noise for the mp — —1 par- 
ticles and vice versa and hence production of pairs from 
our starting 'vacuum' provides a thermalization mecha- 
nism for a given particle species. The density operator 
depends on the wavevector k indicating that the differ- 
ent modes of the system would receive different amounts 
of thermalization, and thus no single temperature may 
be associated to one or the other species for the entire 



system. We may use the above results to produce sim- 
ulations of the system similar to those of Saito et ali^. 
Saito et al. use the classical Gross-Pitaevskii equations 
but build in white noise to represent the quantum fluc- 
tuations. The white noise is chosen to reproduce the 
variance found in the fully quantum equations of motion 
which we use. Here we follow the the Q-function method 
of quantum optics. The Q- function of a density operator 
p is defined as^^ 

Q{a,a*) = - {a\p\a) . (69) 

TT 

For a pure state this is proportional to the probability of 
the state to be in a coherent state a. For our represen- 
tation of state \ f(t)) we may write: 



^|®k((aik|®(a-i-k|)|/(t))|' 



nf^J^exp 

k ^ 



^k ^k I 1 2 I |2 I i+ * * I i+* 

|aik| - |a-i-k| + 0k aik"-i-k + 'Pk "ikO-i-k 



(70) 



(71) 
(72) 



where ij is the number of modes and we have used the 
expansion of coherent states in terms of number states. 



I a ) = exp 



2 



E 



a 
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(73) 



Given this probability distribution we notice that each 
Qk is an exponential of a quadratic form and is thus gaus- 
sian in nature. A change of basis to diagonal form allows 
us to calculate the variance and select a random distri- 
bution of coherent state amplitudes consistent with this 
probability distribution. The random realizations thus 
obtained support our picture of quantum noise providing 
the seeds for the formation of domains, and reproduce 
the results of Saito et ali^. 

The thermal nature of the fluctuations may be verifi- 
able experimentally An experiment could be undertaken 
that considers two settings: ('i^ a system with a starting 
state of N particles that have Bose-condensed in a ground 



state of the system, say with mp^ = 1 (any uniform ro- 
tation of such a state would do, since we are selecting 
one of the degenerate set of ground states for the sys- 
tem). For such an initial state, usually approximated as 
a coherent state, the number variance (or a correspond- 
ing density-density correlato r^^i'^^i^" ) for the particles in 
the mp^ = +1 (or —1) state is Poissonian in nature, i.e. 
CT/ik=o ~ N . The density-density correlator for non-zero 
wavevector q, is calculated for particles in the = 1 
state and corresponds to having Poissonian shot noise: 



a 



+q 



N 



-l{N\pqP-q\ 



N' 



F,=+l 



(74) 



where Pq = '^ik-q'^ik the time independent density 
operator for mp — +1 particles and N'z is the average 
number of these particles in this starting state. The cor- 
relator is independent of wavevector q, as well as being 
time-independent - since we begin in a ground state that 
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Correlation ratio A/, 



Production to source correlation ratio 



100 





y/2rn\cjn 




y/2m\c2\n 



FIG. 9: (Color online) The graphs show the correlation Aq 
for: (a) a range of times and wavevectors and (b) for a given 
time of t — 4 X 16ms. At this time the signal to Poissonian 
noise ratio is 0(10^). 



has no dynamics. 

(a) In the second proposed setting, the starting 
state would be that corresponding to the Berkeley 
experiments^, and modeled in this paper, with N parti- 
cles in the rnp^ — state. This we predict would exhibit 
a thermal distribution for the particles in the mp — +1 
state, showing super-Poissonian statistics. The density- 
density correlator for non-zero wavevector q for this state 
is time dependent, and we use the time dependent cre- 
ation operators already derived in this paper: 



= ^lVlk(l+lVlk-q), 



(75) 



where pq(i) = J2 



,{t)aik{t) is the time dependent 



density operator for mp, — -fl particles, and A'^ik was 
defined in Eg. ([55]). The comparison of the two settings, 
by looking at the ratio of the experimentally computed 
correlations for each case and comparing to the theoreti- 
cal predictions, would provide a direct verification of the 
unique statistics associated with the two-mode squeezing 
in these F = 1 spinor condensates. We may define a 



FIG. 10: (Color online) The ratio Coq/A^o(^) grows above 
unity, at later times (indicated by the region above the black 
curve) . 



correlation ratio deviation, Aq, implicitly, 



Oq 



a 



= 1 



(76) 



+q 



The correlation ratio deviation shows us how much the 
noise differs from that of a coherent state and corresponds 
to a signal to noise figure for an experiment. For our 
model we may calculate numerically Aq as a function 
of time, using for C+q the number of mp = +1 parti- 
cles present in our sample at a given time, having been 
produced from our initial mp = source. 



Aq = 



EA^lkA^lk-q 



Ik 



(77) 



Our calculations, shown in Fig. ([9]), indicate an ex- 
pected Aq to be 0(10^) for initial times i ~ 4 x 16ms. 
This suggests that the thermal noise should be readily ob- 
servable in a typical experimental setting. Furthermore if 
we look at the ratio Coq/No{t), Fig. fTU]) . which is the ra- 
tio of the density-density correlator for our starting state 
to the Poissonian density-density correlator expected for 
the number of mi? = source particles remaining in the 
sample at any given time, No{t), we see that, say at times 
^ 4 X 16ms, this is greater than unity. This suggests that 
at this time, and onwards, the thermal correlations are 
greater than the Poissonian system shot noise and should 
be experimentally observable. Such an experiment to 
probe the correlations in Bosc-Einstein condensed gases 
has been proposed by Altman et al.^*'. The experimental 
effect of super-Poissonian statistics is reminiscent of the 
Hanbury-Twiss-Brown effects^ ^, and the varying number 
statistics in our case are due to the perfect correlations 
between nip = ±1 pairs in the Berkeley experiments. 
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V. SUMMARY 

In this paper we have discussed the seeding of domains 
in the magnetic system of a Bose-Einstein spinor con- 
densate of ^^Rb atoms, in the spin triplet F = 1. The 
analysis was performed in the context of the experiment 
performed by Sadler et alii. The spinor gas behaves ferro- 
magnetically, and the experiment showed the emergence 
of domains of transverse magnetization evolving from an 
initial polar state in the mp — 0, that becomes dy- 
namically unstable as the initial magnetic field is quickly 
ramped down. 

We have described the dynamical quantum fluctua- 
tions of such a sample that starts as a condensate of 
A'^ atoms in a pure F — 1, nip = in analogy to the 
'two-mode squeezing' of quantum optics. Specifically we 
have modeled the initial mp — condensate as a source 
for the creation of particle pairs in the mp = ±1 states. 
Our considerations for the single-mode Hamiltonian were 
extended to a multi-mode approximation which lead us 
to a representation of the system via an su(l,l) algebra. 
Even though the system as a whole is described by a 
pure state with zero entropy, it turned out that consid- 
ering only one species at a given wavevector at a time we 
obtain super-Poissonian fluctuations that correspond to 
a thermal state. This may be thought of as equivalent to 
considering the reduced density matrix for the mp = 1 
degree of freedom, by tracing out the mp = { — 1,0} de- 
grees of freedom. These quantum fluctuations of the ini- 
tial dynamics of the system provide the seeds for the 
formation of domains of ferromagnetically aligned spins. 
The super-Poissonian fluctuations should be observable 
experimentally by looking at density-density correlations 
to compare our time-evolved states to coherent states. 



We should finally note that if we were to use a mean- 
field Gross-Pitaevskii approach then we would not ob- 
serve the evolution of the system in the absence of 
noise. Our approach provides the source of the insta- 
bility needed for the Gross-Pitaevskii equation, in terms 
of the quantum noise of squeezed quantum states. As an 
improvement on the standard mean field description of 
the system, one may try to match the fluctuations of the 
quantum states to the initial boundary conditions for the 
mean-fleld Gross-Pitaevski equations - such an approach 
has appeared in the work of Saito et ali^. This provides 
a semi-classical but fully non-linear description that may 
probe the system in the metastable time regime, after the 
domain formation is complete, to investigate their stabil- 
ity and evolution. The evolution of the uniform start- 
ing state into domains may be viewed a process akin 
to the Kibble- Zurek mechanism^iiS^dSid^ originally de- 
veloped to describe the initial evolution of cosmological 
defects in the early Universe. Furthermore the processes 
described have analogues in different models for the pro- 
duction and statistics of particles in cosmology22il£iS. It 
may be possible in the future to even make the connection 
more explicit and perhaps even get to the point of using 
Bose-Einstein spinor condensates as analogue models for 
cosmological considerations^^'^'^. 
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